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O-t. Abstract 

Some properties of the rho vector meson are calculated within the Nambu- 
Jona-Lasinio model, including processes that go beyond the random phase 
approximation. To classify the higher order contributions, we adopt 1/N C as 
expansion parameter. In particular, we evaluate the leading order contribu- 
tions to the p — > 7T7T decay width, obtaining the value T = 118 MeV, and to the 
shift of the rho mass which turns out to be lowered by 64 MeV with respect 
to its RPA value. A set of model parameters is determined accordingly. 
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I. INTRODUCTION 



Vector mesons are known to play an important role in different aspects of low-energy 
hadronic physics. In fact, effective lagrangians including nucleons, pions and vector mesons 
as fundamental degrees of freedom have been successful in describing a large amount of 
empirical data |IJ. On the other hand, with the introduction of QCD as the fundamental 
theory of strong interactions, it became clear that the structure and properties of both 
baryons and mesons should be understood in terms of quark and gluon degrees of freedom. 
Although much effort has been made in trying to predict low energy hadron observables 
directly from QCD, one is still far from reaching this goal. In such a situation it proves 
convenient to turn to the study of effective models. 

One of the models that has received considerable attention in recent years is the Nambu 
and Jona-Lasinio (NJL) model 0. It was originally proposed as a way to explain the origin 
of the nucleon mass together with the existence of the (almost) massless pion in terms of 
the spontaneous breakdown of chiral symmetry. Later on, as it became clear that quarks 
are the fundamental components of the hadrons, the model was reformulated by replacing 
the nucleon fields by quark fields. In its modern version the NJL model is constructed 

in such a way that all the symmetries (or approximate low-energy symmetries, such as 
chiral symmetry) of QCD are respected. The basic difference with QCD is that gluons are 
eliminated in favor of an effective local quark-quark interaction. For convenient values of the 
coupling constants, chiral symmetry is spontaneously broken, a dynamical, "constituent" 
quark mass is generated and the pion appears as a massless Goldstone boson. The pion 
mass can be brought to its physical value by adding a small bare, "current" quark mass in 
the effective lagrangian. 

Due to the local character of the interaction, the NJL model in 3 + 1 dimensions is not 
renormalizable. Therefore, some kind of regularization scheme has to be introduced. One 
of the methods widely used is to set a cut-off A in the quark 4-momentum. Typical values 
are A « 1 GeV. In a sense, the cut-off gives an idea of the range of applicability of the 
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model. Since the pion mass is well below 1 GeV, one might think that the predictions of the 
model should be reasonable in that sector. On the other hand, the rho mass is not far from 
the typical cut-off values and one might worry about the validity of the model to describe 
the properties of this meson. Within the NJL model, meson modes are viewed as collective 
solutions of a Bethe-Salpeter equation in the random phase approximation (RPA). While the 
pion can be considered as a collective mode, the situation concerning the vector mesons is less 
clear. Since at the mean field level the dynamical quark mass is momentum independent, the 
NJL model does not confine. Therefore, depending on the model parameters, the rho mass 
can be above the unphysical qq threshold and consequently its decay into a qq configuration 
(Landau damping) is possible within the model. In spite of these potential inconveniences, 
at least at the RPA level, the vector mesons seem to emerge as poles in the corresponding 
T-matrix channel with predicted weak decay constants in reasonable agreement with their 
empirical values [[7].|8| . 

To make more definite statements about the validity of the NJL model in the vector 
meson sector, one should study other properties. In particular the width of the rho meson 
in its main decay channel, namely T^^. This requires, however, to go beyond the RPA 
approximation and brings us into the problem of introducing a reasonable classification of 
higher order diagrams. Recently |p|,[lO|| , it has been suggested that a convenient way to do so 
is in terms of an expansion in the inverse of the number of colors, N c . Within QCD, the idea 
of using 1/N C as an expansion parameter has a long history. In Refs. |lT],|l2j it was shown 
that, in order to have a sensible theory for large values of N c , the quark-gluon coupling 
constant has to scale as l/y/N~ c . It follows that, for example, the baryon masses turn out to 
be of 0(N C ), the meson masses of 0(1), and the coupling of a meson to two-meson states 
or to qq states of 0(1/ v/AQ. Independent of these developments, around the same time 
very similar ideas were introduced in a rather different context. Namely, it was shown that 
the traditional Hartree + RPA treatment of many-body problems could be considered as 
the leading order of an expansion in terms of the inverse of the degeneracy of the available 
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Hilbert space, N Moreover, a method to go beyond the leading terms based on the 1/N 
expansion was developed ||14|| . Since the methods applied in the study the NJL model are 
basically the same as those used in many-body physics, it is quite natural to expect that the 
equivalent of the 1/N expansion could be performed. In fact, the results of Ref. JED] seem 
to confirm this expectation. Some previous estimates of the p decay width into two pions 



within the NJL model already exist f!5 |. In general, however, these calculations show that 
this quantity is largely underestimated. As we will see, a consistent treatment in terms of 
1/N C tends to improve this situation. 

In the present paper we report on a calculation of 1/N C corrections to the RPA descrip- 
tion of the p meson within the NJL model. The paper is organized as follows. In Sec. II we 
briefly review the main features of the NJL model. Namely, we show how quark self-energy 
and meson masses, together with their couplings, are obtained in the Hartree + RPA ap- 
proximation and we give the values of the model parameters to be used in our numerical 
calculations. In Sec. Ill we calculate the leading contribution to the p — > nn decay width. In 
Sec. IV we calculate the 1/N C correction to the rho mass. Conclusions are given in Sec. V. 
Some details of the calculations are collected in Appendices A and B. 



II. GENERALIZED NJL MODEL IN THE HARTREE + RPA APPROXIMATION 



In this work we consider the two-flavor version (up and down quarks) of the extended 
NJL model, defined by the effective Lagrangian 

£njl = $(x) (ip - m ) i)(x) + C int . (1) 

In writing Eq.(^J) we have assumed that the bare quark masses are degenerate and set 
m u = ma = mo- The four-fermion interaction term C int is given by 
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Following the usual treatment, we first apply the mean field formalism in the Hartree ap- 
proximation. This leads to the so-called "gap equation" for the quark self-energy (also called 
constituent quark mass) 



m q 
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m + m q GiN c Si 7—4^ j ■ ( 3 ) 

J (2ir) p 2 -m l q 



The integral in Eq.@ is logarithmically divergent. Therefore a cut-off A has to be introduced 
in order to regularize it. In the case of the covariant regularization method, that will be used 
throughout this paper, this integral has an explicit analytical form ||. 

As well-known, when itlq = there is a certain critical value of G\ above which Eq.(0) 
has a non-trivial solution with m q 7^ 0. This corresponds to the dynamical breakdown of 
chiral symmetry. Correspondingly, a non-vanishing value of the quark condensate < ipip > 
develops. For small values of mo, the transition is no longer sharp but the same qualitative 
behaviour is obtained. 

From Eq.(^J), it is seen that the dynamically generated contribution to m q in the Hartree 
approximation scales with GiN c . Therefore, in order to have a sensible large N c limit we 
should impose G\ ~ 0(1/N C ). This type of requirement is completely equivalent to the 
one used in large- A^ c QCD. In fact, the scaling G% ~ 0(1/N C ) was to be expected since the 
effective quark-quark interaction can be thought of as the heavy-gluon (i.e. strongly dressed) 
limit of a gluon mediated interaction. Therefore, in the Hartree approximation m q ~ 0(1). 
It is easy to verify JlO] that the Fock contribution to the self-energy is of 0(1/N C ). 



Making use of the vacuum in the mean field approximation, we shall now study the 
associated meson fluctuations. For that purpose, we solve the Bethe-Salpeter equation for 
the T-matrix 

zT(g 2 ) = iJC - Tr f£* iJC iS(p + iq) x (4) 
J (2tt) 

%T (q 2 ) iS(p - \q) , 

where Sip) = -r— is the fermion Feynman propagator. The T-matrix as well as the kernel 



/C, can be decomposed in terms of scalar, pseudoscalar, vector and axial vector channels. [] 

Due to the presence of the explicit chiral symmetry breaking induced by the current mass 
m , we have to deal with the mixing between pseudoscalar and longitudinal axial fields, called 
7r — A\ mixing. Therefore, in the pseudoscalar channel one obtains a 2 x 2 matrix equation, 
while in the vector channel a 1 x 1 equation is obtained. The solution of these equations can 
be found by summing the corresponding geometrical series. One finally gets 

i%(q 2 ) = (i l5 r a i l5 r a ) — L- x (5) 



where 
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tT p (q 2 ) = ( lfl T a ® lu r a ) x (6) 
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In these equations, j£p are the polarization functions (i.e. quark-antiquark loops) of the 
corresponding channels. They are divergent functions, that require regularization. For con- 
sistency, the same regularization method used in the mean field approximation is to be used 
here. Explicit expressions for these functions can be found in the literature |§. 

The masses of the bound states and resonances are obtained from the position of the poles 
of the T-matrix. In order to determine the meson-gg coupling constants, we approximate 
the T-matrix close to the pole position. We therefore write 



In what follows we will concentrate only on the pseudoscalar and vector channels, since the others 
are not relevant for our purposes. 



i%(q 2 ) « (t l5 T a ®t l5 r a ) x (7) 
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where d = and 
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T p (q 2 )^( lfl T a ® lu T a ) x (8) 



-i{g^ u -q^q u /q 2 ) . 
q z — m z p 

from which the meson-gg vertices can be obtained (see Fig.Q). 

For the vector channel, the rho-quark g p - qq - coupling constant can be defined as the 
corresponding residue at the pole and expressed in terms of the function Jyy. Due to the 
7r — A\ mixing, the situation for the pion is slightly more complicated, but one can still define 
the standard quantities g^qq and in terms of the functions Jpa and Jaa- 

The four parameters of the model can be fixed by fitting four physical quantities, namely 
the quark condensate <qq>( = <uu> = <dd>, since u and d are degenerate), the pion 
and rho masses and the pion decay constant f^. The values we used for the parameters of 
the model are 

A = 1050 MeV GiA 2 = 10.1 

m = 3.33 MeV G 2 A 2 = -14.4 

With this choice, the results in the Hartree + RPA approximation are 

< qq >V3 = _293 MeV m q = 463 MeV 

m w = 139 MeV g^ ql = 4.94 

= 834 MeV = 0.46 

U = 93 MeV g p _ q7j = 2.12 

We point out that the value of the quark condensate, which in the NJL model is simply related 
to the constituent quark mass m q , is consistent with the values obtained in lattice gauge 
calculations. Our prediction for m q , although large with respect to some recent estimates 
HJl6[|, is well within the range of values usually used in the literature |3],f|,[17| • With such 



value of m q the RPA prediction for the p mass, , is below the qq threshold and therefore 
this particle can be clearly identified as a pole in the corresponding T-channel. Note that 
is larger than the experimental value m^ mp = 770 MeV. As shown in section IV, the p 
mass will be lowered to this value by self-energy corrections 2 . 

Before concluding this section it is worth recalling the A^-behaviour of the meson proper- 
ties obtained in the RPA approximation. It can be shown that all the polarization functions 
are proportional to N c . As a consequence, the meson masses are of 0(1), while the meson- 
qq coupling constants are of 0(1/ \fW^}. This is in agreement with the N c counting rules 
obtained in QCD |TT|JT^1 and also with those found in many-body theories [|l3|.|n[ . 



III. THE p -> -kit DECAY 

The main decay channel of the rho meson is into two pions. The leading contribution to 
this process corresponds to the diagram shown in Fig.||, which is of 0{l/y/W c ) in agreement 
with the large N c QCD result. This behaviour results from the factor N c due to the quark 
loop and from the factors l/\fW c associated with the meson-gg vertices (cf. Fig.[l]). 

Using the usual Feynman rules, one obtains the amplitude for the decay process 

M(q;k 1 ,k 2 ) = e tl (q)T^(q;k 1 ,k 2 ) , (9) 

where q = k\ + k 2 , and 

r-fe *,*») = - trfgfr i M x (10) 
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2 Higher order corrections might, of course, also affect the quark and pion properties. However, due 
to the constraints imposed by chiral symmetry in those channels, the calculation of such corrections 
is a quite delicate matter [18|. Here, we prefer to ignore these effects and choose our parameters in 
such a way that the RPA predictions for the pion properties agree well with the empirical values. 



q 

- (Ai <-> fc 2 ) • 



j> + - m g 
Evaluating the trace one gets 

T^q;h,k 2 ) = -2iN c g^ q - q gl_ qq x (11) 
[tf G(q 2 ; kl k 2 2 ) - K G(q 2 ; k 2 , k\)\ , 

where the function G(q 2 ; kf, k 2 ) can be written in terms of the standard regularized functions 
which appear in the NJL model. Its explicit expression is given in Appendix A. 

Making use of Fermi's Golden Rule, the corresponding decay width is found to be 



where fdcfi^ = J ^^ Ek (2^2E k ( 27r ) 4 ^ 4 (ff - fcl _ * s the standard two-body phase-space- 
measure. 

The bar over the squared transition amplitude implies an average over the initial polar- 
izations and a sum over the final ones. Moreover, all particles are put on their mass shells. 
The calculation of expression ([TJ) is straightforward and details are given in Appendix A. 

One obtains 

N 2 ( Am 2 \ 3/2 

G\m 2 -ml). 

In order to be fully consistent with the 1/N C expansion, all the quantities appearing in 
Eq.(p!3|) have to be taken at the values obtained in Hartree + RPA approximation, namely 
those given at the end of the previous section. In particular, we have to take RPA value 
for the p mass. This is equivalent to use the Rayleigh-Schroringer (RS) perturbation theory 
19| . The final result is 



= 118 MeV. (14) 



Q 



As we see, already at leading order in 1/N C we obtain a value which compares reasonably 
well with empirical value [p0| 1 

= 151.2 ± 1.2 MeV. (15) 

It is interesting to note that, instead of using the RPA value, we could have evaluated 
the width using the corrected p mass (i.e. its value once 1/N C corrections are included; see 
next section). This corresponds to use the so-called Brillouin- Wigner (BW) perturbation 
theory [|19|]. This type of expansion implies the inclusion of diagrams of order higher than 
1/y/Nc, in which a p meson line appears as an intermediate state. At least in the case of 

since 



many-body systems, this leads to a poor convergence of the perturbative series J2T 
for a given order in 1/N C , such diagrams usually have opposite signs with respect to all the 
other (disregarded) diagrams of the same order. Within this scheme, the resulting width is 

r^(BW) = 96 MeV, (16) 

which is in somewhat worse agreement with the empirical value. 

We conclude this section with a short remark on the evaluation of T p _, n7T using the naive 
low momentum expansion, as sometimes done in the literature. Within this approximation, 
the form factor G 2 (q 2 ;m 2 ) in Eq. (|l3|) is replaced by its value at q 2 = 0. This is based on 
the assumption that this function is only weakly dependent on q and leads, of course, to an 
important simplification in the evaluation of the corresponding decay amplitudes. In order 
to investigate whether this approximation could have been used in the present calculation we 
display in Fig. (3) the momentum dependence of G 2 (q 2 ; m 2 ) as calculated using Eq. (|A8|) . As 
we see, in the range of momenta we are interested in ( q ~ 800 MeV), G 2 (q 2 ; m 2 ) differs from 
its value at zero momentum for about 30 — 40%. This clearly indicates that, at least within 
the present regularization scheme, the use of the low momentum approximation would have 
been inappropiate. 
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IV. SELF-ENERGY CORRECTION TO THE RHO MASS 



The leading self-energy corrections (i.e. beyond RPA) to the rho meson propagator are 
of order 1/N C , as can be seen from the corresponding diagrams displayed in Fig.^. The 
self-energy, which we will call II Atiy (g), can be decomposed into its real and imaginary parts 
which can be calculated using cutting rules and dispersion relations. 

As already mentioned, the NJL model does not confine, i.e. it allows hadrons to decay 



into free quarks. In order to cure this unphysical behaviour we follow Ref. [22] and impose 
confinement "by hand", that is we neglect contributions from cuts of the diagrams across 
quark lines thus taking into account only cuts across meson lines. This is consistent with the 
pole approximation for the meson lines already used in the previous sections. Within such 
an approximation, the diagrams (b) and (c) of Fig. 2 do not contribute, while only the cut 
across the pion lines in diagram (a) will give a contribution to the imaginary part of FP" (q). 
Therefore 

-\U^{q)= f-%iT»(q;p+lq,p-±q) x (17) 
J (2vr) 



i 



(p - ig) 2 - ml 
i 



iT"{q;p- %q,p + \q) x 



(P + §<?) 2 - m l ' 

where iT^ is the p — ittt amplitude of Eq. ( |10D and it represents the quark triangular loops. 
Making use of the cutting rules, one obtains 

1m W( q ) = -I J d<t>® iT»{q) (iT»(q)T , (18) 

where the pion momenta k% and k 2 are set on shell. Using the transversality property of the 
self-energy, that is U^^q) = {g^ u — q^q u /q 2 ) H(q 2 ), it is possible to write (see Appendix B) 

JmS(g 2 ) = - g r p (g), (19) 

where T p (q) has the same expression as in Eq. (|T2|) except that the rho momentum is not on 
shell. 
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The real part of the self-energy can be calculated from the above expression making use 



It should be noticed that a priori there is no reason to use in Eq. (|20|) the same cutoff as in, 
for example, Eq.(^|). However, if one assumes that the cutoff represents the momentum scale 
up to which the low energy effective model can be used, the fact that certain particles can 
be described up to momenta higher than such scale while others not is quite hard to justify. 
Consequently, we choose to define our model by using the same cutoff for all the possible 
types of loops. 

The pole position is now moved in the complex plane with respect to the RPA value, 
the imaginary part being related to the decay width and the real part with the mass of the 
particle. Then, the expression of the rho mass which includes corrections up to order 1/N C 
reads 



Using Eq.(gjj) together with our set of model parameters, we find a shift of —64 MeV from 
the RPA value. The new, "physical" rho mass is therefore 



which justifies a posteriori the choice made at the end of section II. Here again we have used 
RS perturbation theory to evaluate the mass shift. 

The use of BW perturbation theory would correspond to solve the equation 



of the Kramers-Kronig relation, leading to Q 




(20) 




(21) 



m p = 770 MeV, 



(22) 




(23) 



3 Strictly speaking eq. (|20|) , which becomes exact in the limit A 2 
to Tie £(<? 2 ). 



oo, defines an approximation 
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As we have already stressed this method is not completely consistent with the 1/N C ex- 
pansion. Numerically, however, the resulting shift —60 MeV is not so different from that 
obtained by using RS perturbation theory. In a way this behaviour reminds of that found 
in the study of nuclear giant resonances, where the predictions for the resonance width are 
more sensitive than the shift of its centroid to the perturbation method used. 



The values of the rho mass shift reported above agree well with those obtained in ref. [ 23 1 
using an effective mesonic action. On the hand, arguments based on the assumption that 
the Pq — uj mass splitting ( « 12 MeV empirically) is basically due to the rho mass shift (see 
e.g. ref. |24j]) seem to indicate that our value is too large. It is clear that since the origin of 
such splitting is not yet fully understood this type of argument has to be taken with some 
care. In any case, if one relaxes the constraint imposed above on the rho momentum cut-off 
and reduces it by about 15%, a new self-consistent calculation shows that the shift can be 
reduced down to a value consistent with the po — u> mass splitting. This variation has only 
a minor effect on the predicted rho decay width which is lowered by about 10%. 



V. CONCLUSIONS 

In this paper we have evaluated the rho decay width into two pions and the correspond- 
ing mass shift within the Nambu and Jona-Lasinio model. Since the calculation of these 
quantities requires to go beyond the usual Hartree + RPA approximation, we have intro- 
duced the inverse of the number of colors 1/N C as a good expansion parameter. In terms of 
this parameter, the Hartree + RPA approximation corresponds to the inclusion of diagrams 
of 0(1), while leading contributions to the decay width and mass shift are of 0(1/N C ). In 
order to be fully consistent within 1/N C , we have used the values obtained in the Hartree 
+ RPA approximation to compute the 1/N C quantities. This is equivalent to the use of 
Rayleigh-Schrodinger perturbation theory in many-body physics. 

The calculated decay width turns out to be 118 MeV to be compared with the experi- 
mental value 151.2 ± 1.2 MeV. This corresponds to a 20% accuracy, which is what one could 

13 



expect from a leading order calculation. For the mass shift we have obtained —64 MeV, 
which is of the order of 10% of the RPA value. This result, together with those obtained 



in the quark sector ||10|| , seems to indicate a rather good convergence of the perturbative 
series in 1/N C . In a way, this fact also justifies a posteriori our choice of model parame- 
ters, which was based on the assumption that next-to-leading corrections to the mass shift 
were negligible. Of course, to make more definite statements about the convergence of the 
1/N C expansion other properties as well as higher orders in the expansion would have to be 
calculated. 

In the evaluation of the mass shift we have disregarded the diagrams that contain cuts 
across quark lines. Such an approximation, which is consistent with the pole approximation, 



various 



is certainly a rather primitive way to implement confinement. Quite recently p5 | 
extension of the NJL model in which quarks are confined through a some effective confine- 
ment mechanism have been considered. An obvious extension of our study would be the 
evaluation of the rho meson properties in such type of models. In fact, shortly after the 
present paper was submitted for publication, some work along this line has been reported 

& 



We conclude that the Nambu-Jona-Lasinio model provides an economic and at the same 
time reasonably accurate picture of the p meson, in particular of its decay width into two 
pions. This is true, provided that a single cut-off energy of the order of 1 GeV is introduced 
to regularize the different ultraviolet divergences and that cuts in the diagrams are taken 
only across meson lines, to avoid decay into free quarks. With these approximations, all the 
many-body techniques can be used to systematically explore the consequences of the variety 
of couplings of meson among themselves, as well as with fermions. 
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APPENDIX A: THE p - irn AMPLITUDE AND DECAY. 



As seen in Sec. Ill, the p — > txtx decay amplitude can be expressed in terms of 



T^q-MM) = -2iN c9p _ qq g 2 x 



[fcf G(q 2 ; kf, kl) - h$ G(q 2 ; k 2 2 , k\)\ 



Here 



with h n 



G(q 2 ;kl,kj) = G^(q 2 ;kl,k 2 2 ) + ^[/ 2 (fc 2 ) - 
-2G^\q 2 -klkD] +hl[-I 2 (kl) + G^(q 2 ;klk 2 2 ) + 



+ 



(q 2 -kj + k 2 ] 

2 n 2 



8mJ q 



Jvv(q 2 )], 



2m„ 



a n and 



G^(q 2 ;kl,k 2 2 ) = I 2 (kl) + 



{ 



q 2 - 2{kj + k\ 



X 



- ukD + i 2 {kir 



[q 2 -(k 2 1 + k 2 2 )][l 2 (q 2 )- 
+ (3A; 2 -A; 2 ) (M + M) h tf. fc?> ^ } . 



In the last two equations we have used 



/ 2 (g 2 ) = 4* J 



+ 



(2n) 2 [p 2 — m 2 ] \{p — q) 2 — m 2 ] 



h(q 2 ;klk 2 2 )=At J 



A d A p 1 



x 



[p 2 -m2][(p + /c 2 ) 2 -m2] ' 



2N r 



Jvv(q 2 ) = V (9 + 2m,) W) - 2m 2 / 2 (0) 



where q — k± + k 2 . When pions are on shell, one obtains a much simpler expression 



T"(q) = -2iN c9p _ qq gl_ qq G(g 2 ;m 2 ) (h - k 2 f, 



IS 



where now we have 

G{q 2 -ml 



:i - k 



{ 



h(mi) + 



1-K 



x 



q 2 — Ami 

{q 2 - 2ml) (W) ~ UmD) + 2m%{q 2 ] m*) 
K Jvv(q 2 )y 



{Ah 



2(1 -K 



4 m 2 , 



This pion-on-shell expression agrees with the result obtained in Ref. ||16|| . Note, however, 
the different type of regularization used in that reference. 
To calculate the decay width, one has to evaluate 



A4(m 2 p , ml) 

A=0,+ ,- 

which can be easily done by using the relations 

E #\q)4 x) (a) = -(<r-Wtf 

X- 

and 



(A9) 



A=0,+, 



(A10) 



8tt \ 



Ami 



(AH) 



obtaining as final result the expression given in Eq . (|13|) . 



APPENDIX B: THE RHO SELF-ENERGY 

From Eq.flTTD in the text we have, relabelling momenta, 

n^(g) = -i /|| ( ^ 2 4 (2vr) 4 5\q — fci — k 2 ) x (Bl) 

14 -ml 

i 

k\ — ml 

Cutting the diagram across pion lines, it is possible to obtain the imaginary part of the 
self-energy. This consists in re-writing the amplitude with appropriate rules and putting 
pions on shell with the prescription 

1fi 



1 2ir5(k 2 -ml)6(k°) . (B2) 



k 2 — m 2 (2ir 
Using Eq.fll^), we can write 

1m S(g 2 ) = ~ J \tT»(q)\ 2 . (B3) 



Making use of Eq.( |All|) for off shell p momentum, it is possible to evaluate the integral in 



the rhight hand side. One obtains 



2 / H™2\3/2 



I-^^^^^l-^l x (B4) 



and therefore the result in Eq. ([19]) . 
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FIG. 2. The p—irir amplitude used to calculate the decay width. It corresponds to the function 



T^(g;fci,fc 2 ) as in Eq.(£D. 
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FIG. 3. The 4-momentum dependence of the function G(q 2 ,m 2 ), used in Eq. (13). Here the 
deviation from its value at q 2 = is shown. 
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